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ABSTRACT 

Individual giant radio pulses (GRPs) from the Crab pulsar last only a few 
microseconds. However, during that time they rank among the brightest objects 
in the radio sky reaching peak flux densities of up to 1500 Jy even at high radio 
frequencies. Our observations show that GRPs can be found in all phases of ordi- 
nary radio emission including the two high frequency components (HFCs) visible 
only between 5 and 9 GHz (Moffett & Hankins 1996). This leads us to believe 
that there is no difference in the emission mechanism of the main pulse (MP), 
inter pulse (IP) and HFCs. High resolution dynamic spectra from our recent 
observations of giant pulses with the Effelsberg telescope at a center frequency 
of 8.35 GHz show distinct spectral maxima within our observational bandwidth 
of 500 MHz for individual pulses. Their narrow band components appear to be 
brighter at higher frequencies (8.6 GHz) than at lower ones (8.1 GHz). Moreover, 
there is an evidence for spectral evolution within and between those structures. 
High frequency features occur earlier than low frequency ones. Strong plasma 
turbulence might be a feasible mechanism for the creation of the high energy 
densities of ~ 6.7 x 10 4 erg cm -3 and brightness temperatures of ~ 10 31 K. 

Subject headings: pulsars - PSR B0531+21 - PSR J0534+2200 - Crab pulsar - 
giant radio pulses 

1. Observational Techniques 

Observations with the Effelsberg 100-m radio telescope began on 25 November 2003 
and ended on 28 November 2003. We used a secondary focus cooled HEMT receiver with 
a center frequency of 8.35 GHz, providing two circularly polarized IF signals with a system 
temperature of 25 K on both channels. With a sky temperature of 8 K, and a contribution 
of 33 K from the Crab nebula, the effective system temperature was 70 K. Two detection 
systems were used. First, a polarimeter with 1.1 GHz bandwidth detected total power of the 
left-hand and right-hand circularly polarized signals (LHC and RHC respectively) as well 
as cosZ( LHC, RHC) and sinZ( LHC, RHC). These four signals were then recorded by the 
standard Effelsberg Pulsar Observation System (EPOS). With the Crab pulsar's dispersion 
measure of 56.8 pc cm" 3 , we had an un-dedispersed time resolution of t samp i e = 890 /is, while 
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our sampling resolution was fixed at ~ 640 fis. EPOS was therefore set to continuously record 
data blocks containing 20 periods divided into 1020 phase bins for all four signals (Fig. 1). 

The minimum detectable flux per bin was AS min = 0.117 Jy. Considering the dispersion 
pulse broadening, the detection limit for a single GRP of duration r grp = 3 /is would be 
AS m i n t sam pi e /Tg r p = 25 Jy. Selected GRPs are shown on Fig. 2. 

A wide band ultra-high resolution detection system similar to the one of Hankins et al. 
(2003) was used as well. The two IF signals (100-600 MHz) carrying the LHC and RHC 
information were sampled and recorded with a fast digital storage oscilloscope (LeCroy 
LC584AL). 4 x 10 6 samples per channel were recorded at the rate of 2 x 10 9 samples per 
second. The storage scope was used in a single shot mode. We used giant radio pulses to 
trigger data acquisition by the digital storage scope. The RHC signal was detected, amplified 
and low-pass filtered (3 kHz), then passed to a second scope (Tektronics 744). The second 
scope triggered the digital scope whenever pulses stronger than 3-4 times rms (75-100 Jy) 
were detected. The recorded waveforms were then transferred to disk. This required manual 
intervention and led to a dead time of about 100 seconds after each pulse. 



2. EPOS Data 

About 2.4 x 10 6 periods were observed with EPOS. Presented results are only based on 
a third of all observed rotations. This selection was mostly caused by an insufficient quality 
of the data. Because of the Crab pulsar's weak signal at 8.35 GHz and the strong nebula 
background, ordinary single pulses were not observable with Effelsberg at that frequency. 
It takes about 20 min integration to assemble a mean profile at that frequency. The giant 
pulses come in short outbursts, about 5 to 20 minutes in duration, and appear extremely 
prominent during such burst phases. 

In our subsequent analysis, we set a threshold level of 5 rms = 125 Jy on the sum of 
RHC and LHC for the same phase bin to count as a detection of a giant pulse. More than 
1300 giant pulses were detected that way (Fig. 2) and their arrival phases were computed by 
using the TEMPO 1 pulsar timing package. The data were aligned using the current Jodrell 
Bank timing model. They were found to be a perfect match to the time of arrivals (TOAs) 
obtained at Jodrell Bank before and after the Effelsberg observing session. Upper panel of 
Fig. 3 shows these phases. We found that the giant pulses occur at all those phases where 
the radio components of the Crab pulsar emission are observed (Fig. 3, bottom panel). 



http: / / pulsar.princeton.edu/tempo 
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For all giant pulses, i.e. regardless of their phases, the histogram of their peak strengths 
at 8.35 GHz can be described by a power law with a slope ~ —3.34 ± 0.19 (Fig. 4). This 
result is consistent with the results obtained by others (e.g. Lundgren et al. 1995). Moreover 
in Fig. 5 , we show the phase resolved distributions of four components: two high frequency 
components (HFC1 and HFC2), and the main and inter pulse components (MP and IP). 
Because of limited statistics we could make a model fit only for the IP component, where a 
power law index ~ —3.13 ± 0.22 was found, and it is consistent with the value of ~ —2.9 
presented by Cordes et al. (2004). So far, there is no evidence that the distributions for 
HFC1 and HFC2 differ from each other. However, their slopes seems to be steeper than the 
slope for the IP component. 

The results obtained so far by us suggest that the physical conditions in the regions 
responsible for HFC emission might be similar to those in the main pulse (MP) and inter 
pulse (IP) emission regions. We plan to pursue this question further in the beginning of 
2005, when a new observation campaign at two frequencies, 4.85 GHz and 8.35 GHz, of the 
GRPs of the Crab pulsar is scheduled. 

3. High Resolution Detections 

Because of the burst-like character of the giant pulse emission and the long dead-times 
in our data acquisition, only 150 pulses were observable with the high resolution equipment 
(for example Fig. 6). The peak strength varied from our threshold of ~ 75 Jy to a few 
rare events which even exceeded 1000 Jy. We computed dynamic spectra for the LHC 
and RHC signals by successive Fourier transforms and squaring of the data. Baseline and 
sensitivity corrections were also applied through the use of bandpass averages. The pulses 
were typically 2-5 fis wide and certainly wider than the resolution of our incoherent software 
de-dispersion, 0.7 fis. Hankins et al. (2003) showed the strong temporal variabilities within 
the giant pulses, which are possibly unresolved because of the limited bandwidth of the 
observations. Although giant pulses can be received over a wide bandwidth, the individual 
pulses do not have a uniform spectrum, as seen in Fig. 7. The strongest emission was 
predominantly detected in the upper quarter of the accessible bandwidth. The individual 
pulses consist of ~ 100 MHz-wide clusters of narrow bv ~ 2 MHz spectral lines waxing and 
waning with time (Fig. 7, upper). High frequency features appear earlier than low frequency 
features. Furthermore, if two giant pulses occur in rapid succession (Fig. 6, 7) separated by 
only 100 fis, their spectra are similar though not identical. The maximum emission of the 
leading pulse occurs at higher frequencies than that of the trailing pulse (Fig. 6, bottom). 
The separation of spectral maxima also decreases in the trailing pulse. Intrinsic fluctuations 
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of the emission process or - alternatively - scintillation could be the cause. 

4. Characteristics of the Emission Mechanism 

The maximum flux of S vgrp — 1500 Jy is detected on a typical timescale of T grp — 
3 x 10~ 6 s. With the Crab's distance of d — 2 kpc and accounting for 50 % bandwidth coverage 
(e v = 0.5, Av> = 1.1 GHz) a peak luminosity of L grp = S vgrp e v Aud 2 [1 — cos(r grp / P)] = 
1.6 x 10 25 erg s _1 is obtained. The pulse duration limits the size of the emitting volume 
to V grp ~ (c T grp) 3 = 7 x 10 15 cm 3 and a rough estimate of the dissipated energy density 
yields L grp T grp V~ p ~ 7 x 10 4 erg cm -3 . For our observations the brightness temperature 
is estimated to T b = 4.4 x 10 31 K applying the relation T b = 2nS vgrp T grp d 2 / (kp 2 r grp ) in 
accordance with Soglasnov et al. (2004). Such extreme brightness temperatures accompanied 
by strong polarization indicate a highly coherent radio emission process. We find that the 
GRPs occur at the same phases as the ordinary radio emission, perhaps dominating the 
radio profile at high frequencies. Such coincidence suggests that the GRPs are of the same 
origin as the radio emission at x ~ 20 r ns . At these emission heights, the interaction 
of relativistic particles with turbulent plasma wave packets in the pulsar magnetosphere 
represents a possible source process for the creation of GRPs (Hankins et al. 2003). Strong 
plasma turbulence is the highly non-linear dynamic equilibrium between density fluctuations 
and radiative plasma modes. The purely electrostatic two-stream instability parallel to the 
pulsar dipolar field serves as the trigger mechanism. In the pulsar magnetosphere, electron- 
positron pairs are created in events which are sharply localized in space and time. Each event 
creates an individual plasma shell. Different shell generations have different bulk speeds and 
can interpenetrate the magnetosphere on the way outwards (Jaroschek et al. 2004). In such 
colliding plasma shells, the plasma two-stream instability is excited, which evolves to a fully 
saturated equilibrium of turbulent plasma modes. Self-induced, strong density fluctuations 
in the high-energy beam are coupled to density fluctuations in the low-energy background 
population via radiation pressure. Such coupling leads to significant bunching of particles 
and thereby to coherent emission expected at a frequency of 7 2 f pe , where 7 denotes the 
Lorentz factor in the frame of the slower plasma shell and v pe is the plasma frequency of an 
intrinsically relativistic plasma. 

The mechanism outlined above is attributed as collisionless Bremsstrahlung (Weather- 
all & Benford 1991) and is identified in self-consistent 3D particle- in-cell (PIC) simulations 
(Jaroschek et al. 2004). In these simulations, coherent radiation emission stimulated by 
strong Langmuir turbulence is studied in slab geometry for similar conditions (particle den- 
sities ~ corotation densities of 10 9 cm -3 and Lorentz factors of 7 = 6). Turbulent energy 
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densities reach up to ~ 10 5 erg cm -3 with rise times and fluctuation timescales of a few 
nanoseconds. The emission of a fraction of this energy density is sufficient for the observed 
GRPs. The PIC simulations show that the emitted radiation has a Poynting flux with mag- 
nitude, temporal structure and linear polarization in good agreement with our observations. 
The radiation would last until the supply of high energy particles is exhausted. 
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Fig. 1. — Average pulse profiles of PSR B053f+2f at 8.35 GHz. Twenty periods and 
corresponding polarisation characteristics are shown. 
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Fig. 2. — Train of strong giant pulses at 8.35 GHz (EPOS, Av = 1.1 GHz, sum of the right 
and left handed circular polarisation signals). Results obtained from 6.7 hr of observation 
with the 100 m Effelsberg telescope. 
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Fig. 3. — Upper: pre-fit phases of pulse arrivals in the night of November 27-28, 2003. Data 
from a single scan were taken. Phases and -0.38 correspond to the location of the IP and 
MP, respectively. Bottom: recorded strength and phases of 1318 GRPs. Location of the 
radio components: HFC2: 0.05, precursor: 0.2, MP: 0.3, IP: 0.7, HFC1: 0.9. 
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Fig. 4. — Peak strength distribution of giant radio pulses from the Crab pulsar regardless 
of their phases. For the number of occurrences larger or equal 10 (marked with *) the 
distribution can be roughly described by a power law (S/N) a with index a ~ —3.34 ± 0.19. 
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Fig. 5. — Same as in Fig. 4, but for four phase slices, containing HFC1 and HFC2, MP and 
IP, respectively. Upper panel is for HFC1 (dotted line) and HFC2 (solid line). Bottom panel 
shows MP (dotted line) and IP (solid line) components. The giant radio pulses distribution 
for IP can be described by a power law (S/N) a with index a ~ —3.13 ±0.22 (for the number 
of occurrences larger or equal 10; marked with *). 
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PSR B0531+21 8.35 GHz RHC+LHC MPUR EFFELSBERG 2003-1 1 -27 22:42:43.727903 
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Fig. 6. — Giant pulse with double structure at 8.35 GHz (LHC, RHC), peak strength 
~ 116 Jy. Upper: de-dispersed profiles. Bottom: de-dispersed dynamic spectrum. 
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dynamic spectrum of Crab GRP No. 185 8.35 GHz LHC MPIfR Effelsberg 




Fig. 7. — Spectral structure of giant pulse components. Upper: sequence of spectra of one 
GRP component. Bottom: intensity in two spectral channels separated by 13.6 MHz. 



